Building Gaussian Cluster States by Linear Optics 
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The linear optical creation of Gaussian cluster states, a potential resource for universal quantum 
computation, is investigated. We show that for any Gaussian cluster state, the canonical generation 
scheme in terms of QND-type interactions, can be entirely replaced by off-line squeezers and beam 
splitters. Moreover, we find that, in terms of squeezing resources, the canonical states are rather 
wasteful and we propose a systematic way to create cheaper states. As an application, we consider 
Gaussian cluster computation in multiple-rail encoding. This encoding may reduce errors due to 
finite squeezing, even when the extra rails are achieved through off-line squeezing and linear optics. 

PACS numbers: 03.67.Lx, 42.50.Dv, 42.25.Hz 



Introduction. — The cluster-state model for quantum 
computation 1] is a conceptually interesting alternative 
to the more conventional circuit model 2]. Once a suit- 
able multi-party entangled cluster state has been pre- 
pared, universal quantum gates can be effected through 
the cluster via only single-party projective measurements 
and feedforward. Though originally based upon qubits, 
the cluster-state model can be also applied to other 
discrete-variable systems (qudits) as well as to contin- 
uous quantum variables 0. 

Linear optics represents one of the most practi- 
cal approaches to the realization of quantum informa- 
tion protocols, both for discrete- variable (DV) 0] and 
continuous- variable (CV) implementations p. In the DV 
case, efficient entangling gates cannot be achieved with 
single photons and linear optics. Nonetheless, probabilis- 
tic gates can be applied off-line to an entangled multi- 
photon state that eventually serves as a resource for the 
on-line computation Q. A similar approach uses DV 
photonic cluster states, leading to a significant reduction 
in the resource consumption 0,0- However, the gener- 
ation of the optical cluster states remains highly proba- 
bilistic in this case. 

Although up to six-qubit single-photon cluster states 
have been created via postselection using nonlinear and 
linear optics 0, a possible deterministic, uncondi- 
tional realization of optical cluster states would be based 
on continuous variables. Here, the resources are squeezed 
states of light and the Gaussian cluster states may be cre- 
ated via quadratic quantum nondemolition (QND) inter- 
actions ^J. These interactions, however, cannot be real- 
ized through beam splitters alone. Additional "on-line" 
squeezers are needed for every single link of the cluster 
state, again rendering the mechanism for cluster genera- 
tion rather inefficient with current technology. Moreover, 
the squeezing of the resource states will always be finite, 
inevitably resulting in errors in the cluster computation. 
In this paper, we will address both issues: the avoidance 
of on-line squeezing in cluster-state generation and the re- 
duction of finite-squeezing induced errors in cluster-state 
computation. 



We will show that for any Gaussian cluster state, the 
canonical generation scheme in terms of QND-type inter- 
actions, can be entirely replaced by off-line squeezers and 
beam splitters. Moreover, we propose a systematic way 
on how to build alternative cluster-type states from po- 
tentially cheaper squeezing resources than needed for the 
canonical states. In any of these linear-optics schemes, 
the resource states require correspondingly more squeez- 
ing to compensate for the lack of extra squeezing in the 
beam-splitter network that replaces the QND coupling 
of the cluster nodes. Nonetheless, the main features of 
the canonical QND-made clusters can be preserved. As 
an example, we consider Gaussian cluster computation in 
multiple-rail encoding. This encoding may reduce errors 
caused by finite squeezing. We will see that multiple-rail 
encoding with linear-optics-made clusters, though requir- 
ing supposedly more off-line squeezing for the extra links 
of the larger clusters, can still lead to the same error 
reduction as for the canonical cluster states. 

We define cluster-type states as those multi-mode 
Gaussian states for which certain quadrature correlations 
become perfect in the limit of infinite squeezing , 

(pa-^x 6 )^0, aeG. (1) 
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Perfect correlations uniquely define the corresponding 
graph state (for the DV case, see 12J ) . Here we use the 
dimensionless "position" and "momentum" operators, x 
and p, corresponding to the quadratures of an optical 
mode with annihilation operator a = x + ip. The modes 
aeG correspond to the vertices of the graph, while the 
modes b G N a are the nearest neighbors of mode a. 

Canonical cluster states via linear optics. — The canon- 
ical way to build CV cluster states would be to send a 
number of single-mode squeezed states through a corre- 
sponding network of QND gates 0, 0] • Each individual 
QND gate could be realized via two beam splitters and 
a pair of on-line squeezers 0. However, by including 
the initial single-mode squeezers into the QND network, 
the resulting total circuit corresponds to a big quadratic 
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Hamiltonian applied to a number of vacuum modes. This 
total transformation can be decomposed into a linear- 
optics circuit followed by single- mode squeezers and a 
second linear-optics circuit [l3|, where the first linear- 
optics circuit has no effect on the vacuum modes. Even- 
tually one has just one linear circuit applied to a number 
of single-mode squeezed states; in principle, this works 
for any cluster (or graph) state. Let us explain this in a 
little more detail. 

The canonical generation of CV cluster states [ll| 
from momentum-squeezed vacuum modes via QND-type 
interactions can be described by x' a = x a and p' a = 
Pa + Y.beN a where x a = e +r xi 0) and p a = e _r pi 0) , 
Va G G, with vacuum modes labeled by superscript (0). 
Note that, according to Eq. (QJ, the canonical cluster 
states become perfectly correlated for r — > oo. We have 
the following linear Bogoliubov (LUBO) transformation, 
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Using = x\ yJ) +ip\ yJ) , we can extract the LUBO matrix 
elements for canonical cluster generation: A aa = coshr, 
B aa = sinhr, A ah = B ah = (i/2)e+ r (V6 G N a ), and 
A a i = B a i = [V/ N a (l ^ a)]. Now the input 
modes of the LUBO transformation are vacuum modes 
instead of the squeezed input modes of the QND net- 
work. The next step is to decompose this LUBO trans- 
formation into a linear-optics circuit V\ a set of single- 
mode squeezers, and another linear-optics circuit U 
Using the singular value decomposition for A and £?, 
we have to satisfy V = A^UA' 1 = B T U*B~ 1 and 
U = AVA^ = BV*B~ D 1 for the unitary matrices U and 
V. Here A 2 D is the diagonalized version of AA^ and A^A; 
similarly, we use Bp for BB^ and (B^ B) T . The singular- 
value matrix equations lead to the conditions for J7, 

TmUai-Ci(r) ^ ReU bl = , Va,/eG, (3) 
beN a 

with 

ReU al [Di(r) - M Na ] - ^ ^ ReU kl = , (4) 

beN a k£N b ,k^a 
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where 



Ci(r) 



Xf sinh r + J Xf cosh r 



(5) 



Di(r) = (Af cosh 2 r-Xf sinh 2 r)/[e +2r (A^ - Af )/4], and 
ReU a i = Re(U a i) : etc.; Mjv a represents the number of 
nearest neighbors of mode a. The double sum in Eq. (@J) 
sums over all second neighbors of mode a (including mul- 
tiple counting of identical neighbors of the nearest neigh- 
bors of mode a) . The expressions J Xf are the singular 
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FIG. 1: Different generation schemes for a canonical two- 
mode cluster state. The input modes are always vacuum 
modes; a) canonical generation via QND interaction; b) QND 
gate replaced by a beam splitter (BS), two on-line squeez- 
ers, and another beam splitter ^^; c) generation of a stan- 
dard two-mode squeezed state (via two squeezers and a beam 
splitter) plus two local squeezers [3; °0 simplest generation 
scheme with the QND gate replaced by a beam splitter. 



values of A and similarly for B. Equation © can be 
incorporated into the column vectors of U, 
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I au + iCi(r)J2 bGNi a b i \ 
&2i +iQ(r) ^2beN 2 a bi 

\a N i + iCi(r)Y, be N N a bi J 
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using ReUki = otku with the constraints of Eq. Q. and 
J2i(ui)k(ui)l, = 5kk' for unitarity. These conditions au- 
tomatically satisfy A = UA D V ] and B = UB D V T , 
where the diagonal matrices Ajj and Bjj contain the 
corresponding singular values. Thus, we effectively con- 
structed a linear-optics circuit U that exactly outputs 
the canonical cluster states when applied to the off- 
line squeezed input modes with squeezed quadratures 

The 
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are "redefined" vac- 



uum modes after the first linear-optics circuit that 
has no effect. 

Let us consider the example of the canonical two-mode 
cluster state. It corresponds to two momentum-squeezed 
modes (squeezed by r) coupled through a quadratic QND 
gate e 2%x ® x , see Fig^i. In this case, we have Xf = 
cosh 2 r + e +2 74, Af = sinh 2 r + e +2r /4, d(r) = C 2 (r) = 
C _1 (r), and Di(r) = 1. Thus, in the equivalent linear- 
optics scheme (FigJIJl), two equally squeezed modes are 
combined at a beam splitter described by U with column 
vectors as in Eq. ©, Z = 1,2, TV = 2; the constraints in 
Eq. are always satisfied. 

A possible solution for J7, choosing otyi = OL21 — 0, is 



U = 



C{r) i 
i C(r) 



(7) 



Each input mode is momentum-squeezed with R\ = 
In (y^Af + R- For r = 0, we obtain Xf 5/4 

and Af = 1/4, and thus R = ln[(l + v / 5)/2]. This is the 
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residual squeezing coming from the QND gate, which is 
now additionally applied off-line before the beam split- 
ter. Note that the QND gate here is not simply replaced 
by a beam splitter followed by two single-mode squeez- 
ers and another beam splitter all together applied 
to two initial squeezed states (FigHJ)). Instead, the cir- 
cuit is further simplified, and only one pair of squeezers 
followed by one beam splitter operation act upon initial 
vacuum modes. However, in order to produce exactly the 
canonical state, this circuit is rather wasteful in terms 
of squeezing resources. The most economical two- mode 
entangled state is a standard two-mode squeezed state 
(TMSS) built from two single-mode squeezed states with 
a 50/50 beam splitter. Like any pure Gaussian two- mode 
state 14], also the canonical two-mode cluster state can 
be obtained from a TMSS via local Gaussian transfor- 
mations, including local squeezers (Fig (HO- 

Comparing resources, for example, an excess noise of 
one vacuum unit in the quadrature correlations p\ — X2 
and p2 — x\ can be achieved with a TMSS built from 
two 3 dB squeezed states (see below). The canonical 
two- mode cluster state with these correlations (r = 0), 
would require two 4.18 dB squeezed states combined at 
the asymmetric beam splitter in Eq. Q. In general, 
the canonical TV-mode cluster states are always biased 
in x and p 0, (x 2 a ) = e +2r ^ (p 2 a ) = e~ 2r + M Na e+ 2r , 
unnecessarily requiring extra local squeezing to achieve 
a certain degree of correlations and entanglement; they 
cannot be obtained from iV-mode pure Gaussian states 
in standard form 16] without the use of local squeez- 
ers. In the following, we shall consider a whole family of 
states producible via linear optics and exhibiting cluster- 
type correlations. This family will include states cheaper 
than the canonical cluster states. 

Cluster-type states via linear optics. — In order to cre- 
ate states with correlations as in Eq. via linear optics, 
we consider p-squeezed input modes, hi 



+Ri -(0) 



ie 



-Rid*) 



p\ , and a general linear-optics transformation, 
a' k = UkiCLi, with a unitary matrix U. Using a' k = 
x' k + ip' k , one obtains the output quadrature operators, 
and in order to satisfy the correlations in Eq. (QJ, we as- 
sume the large noise terms (those proportional to e +Rl ) 
cancel. This is possible if and only if 



times one unit of vacuum quadrature noise (1/4 in our 
scales) which we omit in the following. These excess 
noises will lead to the errors in the cluster computation 
Q. Let us consider two- mode cluster- type states. 
Look at the simple circuit described by 
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applied to two momentum-squeezed modes. The circuit 
is a 50/50 beam splitter with Fourier transforms of mode 
2 before and after the beam splitter. It clearly satisfies 
the conditions in Eq. JBJ). The resulting state is the stan- 
dard TMSS up to a local Fourier transform of mode 2. 
No extra local squeezers are needed to obtain this cluster- 
type state from standard two- mode entanglement. Ac- 
cording to Eq. (0), the excess noises in the cluster- type 
correlations are 2e~ 2Rl in pi — x 2 and 2e~ 2R<2 in p 2 — £i, 
corresponding to the optimal entanglement for a given 
input squeezing. Note that the solutions in Eq. (fTUj) and 
Eq. coincide for r — > oo when C(r) — > 1. 

Incorporating Eq. (JSJ into the column vectors of U 
leads to 



ui = 



/ an + i ^2 hGNl a bi \ 
\a m +i Y,beN N a bi J 



(11) 



with J2i(ui)k(ui) k , = Skk' fc> r unitarity and again 
KeUki = otki- Using this formalism, for any given graph 
state, one obtains a set of geometrical conditions for N 
real vectors a k = (a^i, 0^2, «/cAr). For instance, a 
linear 4- mode cluster state has aids = 0*20*4 = —1/5, 

0102 = d\d^ — 02 03 = 0304 = 0, I \di 1 1 = 1 1 ^ 4 1 1 = 3/5, 

and \\02W = 1 1«3 1 1 = 2/5, achievable, for example, 
with ol = (1/^,1/^10,0,0), 0% = (0,0,-2/^10,0), 
ol = (0,-2/^,0,0), andaj = (0, 0, 1/VT0, -1/v^). 
This, together with Eq. ([TTj). where we set N = 4, 
defines a possible solution U. In general, any 4-mode 
transformation U can be decomposed into a network of 
4(4 — l)/2 = 6 beam splitters [Tjj- However, for the 
linear 4-mode cluster, there is a minimal decomposition 
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After inserting these conditions into the excess noise 
terms (those that vanish for infinite squeezing R\ — » 00), 
we find that every quadrature correlation of Eq. (QJ la- 
beled by a has an excess noise variance 



E 



ReU al (1 + M Na ) + J2 £ Ret/, 



beN a keN b ,k^a 



-2Ri 
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FIG. 2: Examples of cluster- type states producible via linear 
optics; each optical mode/cluster node a is represented by a 
vector a a (see text): a) linear 4-mode cluster, nonlinear 6- 
mode cluster; b) "EPR/GHZ-type" clusters: linear 2 and 
3-mode clusters, T-shape and cross-shape clusters. 
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a) c) 




FIG. 3: Error propagation when an input state is attached 
to and teleported through a) a linear 3-mode cluster; b) a 4- 
mode diamond cluster; c) an iV-mode multiple-rail cluster; all 
clusters are generated via off-line squeezing and linear optics. 



of U into a network of only three beam splitters U = 
F 4 S 12 FiB+ 4 (l/V2)B? 2 (-l/V2)Bz 3 (l/V5)F^Fl Here, 
F k (Fourier transform of mode k) is the 4-mode iden- 
tity matrix h except (F k ) kk = h B^(t) (beam split- 
ter transformation of modes k and I) equals I4 except 
(Bp kk = t, (B±) kl = VT^, (B±) lk = ±VT^V, and 
(B kl )u = =pt, and the matrix S\ 2 swaps modes 1 and 2. 

Similarly, one can find linear-optics solutions for other 
cluster- type states (FigE). The example of a nonlin- 
ear 6-mode cluster in FigEK ( a potential resource for 
two-mode evolutions) leads to two orthogonal subspaces 
{a/c,d^,a m } with {fc,Z,ra} = {1,3,5} and {k,l,m} = 
{6,4,2}, where a k a t -3/10, d k d m did m -1/10, 
and 11^11 = 11^11 = 7/10, ||a m || = 3/10. A possible so- 
lution is (1/V5, 1/V2, 0, 0, 0, 0), (l/vt, -1/V2, 0, 0, 0, 0), 

(-1/^20,0,1/2,0,0,0) for vectors {1,3,5}, respectively, 
and also (0, 0, 0, -1/V^O, 0, 1/2), (0, 0, 0, 1/VE, l/y/2, 0), 
(0,0,0,1/^5,-1/^2,0) for {2,4,6}. In general, using 
the geometrical conditions for N TV-dimensional vectors, 
solutions can be recursively constructed for any graph 
state: start with aj = (an, 0,0, . . .) determined by the 
norm ||a?i||, then = (a^i, a 22 , 0, . . .) given by the over- 
lap d\d 2 and the norm ||d?2||, dj = (#31, 0^32, <^33, . . .) via 
d\ds, d 2 ds^ and ||c?3||, etc. 

Redundant encoding for error filtration. — As an appli- 
cation, we consider multiple-rail encoding in a Gaussian 
cluster computation protocol |2|. Look at the diamond 
cluster state in Fig|2t>. In this case, the vector con- 
ditions become d\ a 2 = aids = d 2 d± = as a 4 = 0, 
a 2 ds = 6l\6l± — —2/5, and ||afc|| = 3/5. This can be 
satisfied, for example, via d[ = 

dl = (0,0,l/V3,-2/\/l5), dl = (0,0,0,^3/5), and 
dj = (0, \/3/5, 0, 0). A simple cluster protocol would 
be to teleport an input mode onto mode 1 of the di- 
amond state and from there to mode 4 (FigUJb). For 
simplicity, we assume that the input mode is attached 
to mode 1 through the QND gate e 2tXin ® Xl such that 
Pi — > Pi + ^in and 

Pin Pin H~ ^1, %"l 

and fin re- 
main unchanged. Now p-measurements of the input 
mode and modes 1 through 3 (with results Si n , S1...3) 
would transform mode 4 into the input state up to some 
known corrections and excess noise coming from the im- 



perfect diamond state. After the correction operation, 
F t X(-5 in )F t X(-5i)F t X(-(52 + 5 3 )/2), with an inverse 
Fourier transform operator acting as x — > p and p — > 
—x, and ^-displacements X(s) such that x — > x + s, mode 
4 is described by x out = X[ n + (pi — x 2 - x 3 ) - (£4 — x 2 — x 3 ) 
and p ou t = pi n -(p 2 -xi-X4)/2-(p3-x 1 -X4 : )/2. Here, 
Pi n and pi correspond to the momentum operators of the 
input mode and mode 1 before the QND coupling. 

We see that the imperfect teleportation fidelities de- 
pend on the correlations of the diamond cluster, Eq. (QJ 
with FigEb. In the corresponding excess noise variances, 
Eq. (0, let us consider only the effect of finite R3 and 
i?4, assuming R\ ^> 1 and R 2 ^> 1. This leads to neg- 
ligible excess noises, a = 1 and a = 4, in £ ut, because 
the vectors d\ and live in the two-dimensional sub- 
space I = 1,2. In Pout, according to Eq. ® with a = 2 
and a = 3, we find a total excess noise of (13e~ 2Rs + 
5e -2i? 4 )/ 12 = 3e~ 2R /2, assuming R 3 = R 4 = R, which 
is half the excess noise compared to any linear 3-mode 
cluster protocol (FigEJt) with two modes highly squeezed 
and one finitely squeezed by R; the generality of this re- 
sult can be easily proven using Eq. (0 and Eq. (JTTJl for 
the 3-mode protocol, with p out — p- XYL — (p 2 — x 1 — x 3) and 
^out = ^in + (jpi-x 2 )- (p3 ~ X2) ~ #in- More generally, 
we obtain a p excess noise of 3e~ 2R /m for multiple-rail 
encoding (FigEt), where m is the number of intermediate 
nodes (number of rails) between mode 1 and the output 
mode As a result, the excess noise in cluster compu- 
tation may be reduced by teleporting the input through 
multiple paths. Remarkably, the linear-optics- made clus- 
ters achieve the same error reduction as obtainable for 
QND-made clusters with QND gates freely available. 

Conclusion. — We showed that any Gaussian cluster 
state can be built via off-line squeezing and linear optics 
without QND couplings. Simple vector conditions lead to 
potentially cheaper cluster- type states. In multiple-rail- 
encoded cluster computation, the same error reduction as 
for the QND-made clusters can be achieved. Our results 
pave the way for experimental realizations of small-scale 
cluster computation with continuous variables. 
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